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Abstract 

As airplane designs have trended toward the expansion 
of flight envelopes into the high angle of attack and 
high angular rate regimes, concerns regarding modeling 
the complex unsteady aerodynamics for simulation 
have arisen. Most current modeling methods still rely 
on traditional body axis damping coefficients that are 
measured using techniques which were intended for 
relatively benign flight conditions. This paper presents 
recent wind tunnel results obtained during large- 
amplitude pitch, roll and yaw testing of several fighter 
airplane configurations. A review of the similitude 
requirements for applying sub-scale test results to full- 
scale conditions is presented. Data is then shown to be 
a strong function of Strouhal number - both the 
traditional damping terms, but also the associated static 
stability terms. Additionally, large effects of sideslip 
are seen in the damping parameter that should be 
included in simulation math models. Finally, an 
example of the inclusion of frequency effects on the 
data in a simulation is shown. 


Nomenclature 

b wing span 

c mean aerodynamic chord 

Q rolling moment coefficient 


Cl f 

ec- 


static lateral stability derivative 
body axis roll damping derivative 


Cj rolling moment due to yaw rate derivative 


stability axis roll damping derivative 
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pitching moment coefficient 
pitching moment at a = 0 
static pitching moment derivative 
pitch damping derivative 
yawing moment coefficient 

Strouhal number (reduced frequency) col/2V 

characteristic length (b/2 or c/2) 

body axis roll rate 

nondimensional roll rate, pb/2V 

body axis pitch rate 

nondimensional pitch rate, qc/2V 

free stream velocity 

angle of attack 

mean angle of oscillation 

angle of sideslip 
roll angle 

oscillation frequency 
stability axis rotation rate 


Introduction 

At high angles of attack, unsteady aerodynamic effects 
may have a major impact on the maneuverability and 
controllability of an airplane. Currently, some modern 
fighter aiqfianes are capable of performing transient 
maneuvers involving high pitch rates to extreme angles 
of attack; this has been graphically demonstrated in the 
so-called "Cobra maneuver" flown by Soviet aircraft at 
recent airshows. The advent of innovative high-a 
control effectors such as thrust vectoring and forebody 
controls will enable even greater capability to 
effectively exploit a substantially enlarged envelope for 
air combat. The impact of unsteady aerodynamic 
effects on airplane flight dynamics (including stability 
and control effects during rapid high angle-of-attack 
maneuvers) and a practical means of utilizing these 
unsteady effects need to be addressed. 


Because current generation developmental aircraft 
programs are relying more heavily on simulation based 
predictions of airplane dynamics to reduce risk and cost 
of flight test programs, accurate prediction of aircraft 
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dynamics is critical. With airplanes operating in the 
highly nonlinear flow regimes, including post stall 
angles of attack, with substantial angular rates, 
prediction of departures and flight dynamics have 
received increased attention. Several efforts have been 
undertaken to address this area, including guidelines 
development for preliminary design 1 , improved testing 
techniques 2 , and improved analysis techniques such as 
in the prediction of falling leaf motions 3 . Nevertheless, 
no industry standard or accepted procedures exist which 
provide for high confidence prediction of airplane 
dynamics in the high angle of attack, and high angular 
rate environment. Some efforts have been made which 
seem to improve the simulation predictive capabilities 4 , 
however these techniques still do not always provide 
good simulation predictions. Most simulation 
development programs today rely on static wind tunnel 
measurements, rotary balance measurements and forced 
oscillation measurements to develop the aerodynamic 
model. 

A significant amount of research has been conducted in 
the area of unsteady aerodynamics at high angles of 
attack. Results for two-dimensional airfoils show 
substantial force overshoots for pitching and plunging 
airfoils 3 ’ 7 . Wind tunnel experiments have been 
conducted using delta wings undergoing pitching and 
plunging oscillations to determine vortical flow 
characteristics during large-amplitude motions 8 '" 1 . The 
primary focus of the studies was the dynamics of the 
leading-edge vortex while the wing was undergoing 
motion. Under these conditions, the vortex burst point 
location was observed to lag the static location during 
pitch-up and pitch-down motions. Water tunnel 
results 11 have shown vortex lag times of up to 30 
convective time units (time required for an air particle 
to travel across the wing), which is a much slower 
response than that normally associated with two- 
dimensional dynamic stall phenomena. Research was 
conducted on a series of wings of different sweep 
angles and showed that the magnitude of the unsteady 
aerodynamic effects decreased as the sweep angle 
increased 12 . Additional tests of semi-span models with 
various wing sweeps were reported in reference 13. 
Reference 14 extended this work to a representative 
three-dimensional configuration, tested at zero sideslip 
and with no control deflections, which showed 
significant force and moment increments due to 
dynamic stall effects. Measured persistence times of 
the dynamic effects were also presented. Reference 1 5 
presented detailed wing surface pressure measurements 
of a straked-wing fighter model undergoing pitch 
oscillations. Force and moment measurements from 
water tunnel studies with a series of straked-wing 
models have also been reported 16 . Recent studies have 
been reevaluating data reduction techniques from 


conventional wind tunnel tests 17 , and mathematical 
structures for simulation modeling 18 " 23 . 

This paper presents results from tests of an F-16XL 
model in the NASA Langley 12 Foot Low Speed Wind 
Tunnel undergoing pitching motions, and another F- 
1 6XL along with a modern fighter model in the 
Langley 14 -by -22 Foot Subsonic Wind Tunnel 
undergoing roll and yaw oscillations at various angular 
rates and amplitudes. These results can be used to help 
develop better understanding of the complex dynamic 
aerodynamics and in turn aid in the development of 
better mathematical modeling techniques. 

Models and Test Apparatus 

Pitch axis static and dynamic tests were conducted in 
the NASA Langley 12-Foot Low- Speed Wind Tunnel 
using a 0.10-scale F-16XL model (figure 1). The 
model was mounted on a dynamic test rig through a six- 
component strain-gauge balance. The dynamic test rig 
is a computer controlled, hydraulically actuated system 
that was sting-mounted on a C-strut support system. 

The mounting arrangement rotated the model about the 
reference center of gravity location of 0.558c, over an 
angle of attack range of -5° to 80°. The maximum 
capability of the dynamic test rig was 260 deg/sec pitch 

rate, and 2290 deg/sec^ pitch acceleration. Further 
description of the dynamic test rig may be found in 
reference 2. The tests were conducted at a dynamic 
pressure of 4 psf resulting in a Reynolds number of 10^ 
based on the mean aerodynamic chord. 

All data were obtained with the leading-edge flaps at 0° 
deflection. Static data were obtained for angles of 
attack from -5° to 80° at zero sideslip and zero 
deflection of trailing-edge surfaces (flaperons and 
elevons). Additional data were obtained for four values 
of sideslip, (p = -5°, -10°, -20°, and -30°), and for two 
values of symmetric trailing-edge control surface 
deflections (8 e = -20° and 20°). Oscillatory data were 

obtained at initial angles of attack (a 0 ) between 20° and 
60° at several frequencies and amplitudes. The effects 
of sideslip and control surface deflections were 
measured only at an initial angle of attack of 35° and 
two frequencies (0.6 and 1.41 Hz) at the same sideslip 
values and control deflections as in the static test. 

Ramp tests (constant pitch rate motions) were 
conducted similarly with various positive and negative 
pitch rates, and various starting and ending angles of 
attack. 

Data recorded during the test runs included a linear 
variable differential transformer reading from the 
dynamic test rig to determine the pitch angle, six- 
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component force and moment data from a strain-gauge 
balance, and wind tunnel dynamic pressure. Data were 
sampled at 1 00 Hz with an in-line 1 00 Hz low pass 
filter. All data channels were subsequently filtered 
using a 6 Hz low pass filter. The combination of model 
pitch rates and tunnel speeds allowed for a realistic 
representation of full-scale maneuvering conditions 
within the capabilities of current and projected future 
high performance airplanes. 

Roll and yaw axis forced oscillation data were obtained 
with an 0.18-scale F-16XL model and other fighter 
models using a forced oscillation rig in the 14 -by -22 
Foot Subsonic Wind Tunnel (figure 2). Data were 
obtained for 4 amplitudes of oscillation and a range of 
frequencies to cover a broad range of angular rates and 
reduced frequency combinations. Data were sampled at 
200 Hz with an inline 1 00 Hz anti-aliasing filter and 
further filtered during post processing with a low pass 
filter similar to the pitch data discussed above. 
Ensemble averaging was conducted over approximately 
40 cycles of the data before using it to develop 
regression models. Time history format of the data 
were obtained and analyzed in addition to in-phase and 
out-of-phase coefficient data used in traditional forced 
oscillation testing. 

Similitude Considerations 


documented in the 1950’s, however this effect typically 
is not addressed in current testing or simulation 
modeling. The only test method currently in use to 
obtain k effects is the forced oscillation test, however 
this data typically is modeled as a linear derivative in 
non-dimensional rate for only one k value which is 
usually picked based on estimates of Dutch-roll 
frequency or short period frequency. Although this 
method worked well for airplanes in the past, the new 
realities of operating in very non-linear flow regimes 
may have invalidated some of the assumptions and 
simplifications used previously. 

A rigorous modeling approach might look something 
like: 


q = f(a, p, 8 , M, Re, Ql/V, k) 


where each coefficient is dependent on all of the 
nondimensional parameters. The current modeling 
approach looks more like: 


Cj = 


f (a,p,8) (measured static data) 

+ f (M) (Mach number effects) 

+ f (Re ) (Reynolds number corrections) 
+ f (Ql/V) (rotary balance data) 

+ f (pl/V, ql/V, rl/V) (forced oscillation data) 


Full-scale predictions using sub-scale models rely on 
satisfying a set of similitude requirements. A summary 
of some of these requirements is shown in figure 3. As 
can be seen from the figure, accurate predictions of full- 
scale characteristics require correlation of many, 
sometimes conflicting, similitude parameters. When 
predicting static characteristics, for example, flow 
angles and control surface deflections need to be 
identical between sub-scale and flight. Mach number 
and Reynolds numbers also need to be maintained 
constant between model and flight for the data to be 
valid as tested. When the aiqtlane is undergoing 
rotational motions, reduced angular velocity and 
Strouhal number (frequency content) need to match 
between model and flight according to the 
nondimensional analysis results. Finally, for 
unconstrained model testing, in order to achieve the 
same (scaleable) flight dynamics, the model needs to be 
dynamically scaled as shown. 

Various types of wind tunnels and test techniques have 
been developed largely to address the many similitude 
requirements. Traditional wind tunnel testing has tried 
to take into account all of the static model 
considerations, and also the angular velocity scaling. 
The frequency-dependence in the data denoted by the 
Strouhal number (k, or reduced frequency) was 


where effects are assumed independent, and also that 
the linear superposition principles apply. This 
approach, although in standard use, is based on broad 
assumptions, which at least until recently, appear to 
have been justified. 

Pitch Axis Results 

Several approaches have recently been applied to 
improve modeling the observed nonlinear, unsteady 
aerodynamic responses seen during large amplitude 
pitching motions. To date, there is no consensus of 
approaches to model observed effects. Little validation 
of proposed modeling methods is available, and the 
literature is dominated by theoretical and mathematical 
studies. Figure 4 shows some basic characteristics 
observed in oscillatory motions. The data show the 
effect of oscillation amplitude on lift and pitching 
moment for constant maximum pitch rates. As can be 
seen, large variation of aerodynamic parameters is 
apparent due to the amplitude. At a = 36°, (midpoint 
of oscillation cycle), all normal state variables (a, p, q) 
are identical; however, the aerodynamic coefficients are 
very different. 

Results from constant pitch rate motions are shown in 
figure 5. These data were all obtained from ramp 
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motions beginning at a = 0° (positive pitch rates) or 
a = 75° (negative pitch rates). The coefficient 
increments were obtained by subtracting the quasi-static 
results (qc/2V = 0.0004) from the dynamic results. 
These increments show large variations in lift and 
pitching moment which are nonlinear functions of a 
and non-dimensional rate. The largest increments seen 
during the pitch-up motions were at the highest angles 
of attack. The largest increments seen (especially in 
lift) during the pitch-down maneuvers were for the 
moderate angle of attack range near the steady stall 
angle of attack. This phenomenon is consistent with 
flow lag effects seen and reported in numerous previous 
studies. At positive pitch rates the flowfield formation 
dynamically lags the static conditions for a given angle 
of attack and results in a delay of the separated flow 
fields. This gives rise to the large increments in lift at 
the high angles of attack and pitch rate. At negative 
pitch rates, again the flow formation lags, such that the 
re-establishment of vortex systems and attached flow 
occurs at lower angles of attack than would be the case 
for static conditions. 

The effect of oscillation amplitude is shown in figure 6. 
The data show that there are large effects of amplitude 
at the negative pitch rates, which increase with 
amplitude. Positive pitch rate results are relatively 
minor with amplitude changes. This illustrates the 
importance of the “time history effect" which is why 
modeling dynamic data without either time dependence 
or new state variables may lead to erroneous results. At 
any given vertical crosscut on figure 6, the standard 
state variables (a, q, & ) are equivalent; however, large 
discrepancies in the data exist. The difference is that 
for each successive increase in amplitude, the model 
previously was further into vortex-dominated flow or 
separated flow regimes before arriving at the given 
angle of attack. The flow lags mentioned previously 
then contribute to the differences seen in the data. For 
positive pitch rates through a = 2 1 °, little change is 
seen with amplitude. This is because the model is 
arriving at the a = 21° point from lower, primarily 
attached flow conditions. The time scales associated 
with flow lags during attached flow are very fast, and 
do not result in significant increments. Figure 7 shows 
similar data at a = 51°. At this angle of attack, the 
large increments are observed for the positive pitch 
rates. Again, similar to the previous example, the small 
amplitude shows the smallest increments. At a = 5 1 °, 
and an oscillation of ± 5 deg, the model remains in 
separated flow with the leading-edge vortex system 
burst. Oscillations inside a completely separated flow 
regime without leading-edge vortex effects result in 
very little dynamic increments. When the amplitude 
increases, and vortex formation occurs during the 


oscillation cycle, then the dynamic effects become 
significant and large increments are seen. This 
behavior is very similar to that discussed in previous 
reports regarding crossing of bifurcation points. 

The experimental study did not conduct sufficient 
testing to make extensive comparisons of amplitude and 
rate effects with constant Strouhal numbers (for 
similitude requirements), but linear derivative estimates 
could be calculated for pitch damping and then 
compared with ramp data for trends. A postulated 
mathematical model; 

C m (a,k) = C mQ (k) + C ma (k)*a + C mq (k)*qc/2V 

was used and considered to be a function of k. 

Parameter estimation using the method of reference 24 
was conducted using this model. The fit of the data, for 
constant k values was very good and generally 
accounted for over 98 percent of the variation seen in 
the data. Figure 8 shows a comparison of the dynamic 
increments, computed as previously discussed, for 
constant rate ramp data with the results obtained from 
the regression model for the C m (k) term at a = 36°. 

The data show a large dependence on k. At large 
values of k, stable damping characteristics are seen. As 
k decreases, the pitch damping stability decreases and 
becomes unstable at the low values of k. The trend 
appears to be consistent with the ramp results which 
show unstable damping at this angle of attack. 

Limited testing at sideslip was also conducted during 
the large amplitude pitching tests. Figure 9 shows the 
general character of results in sideslip. The lags in the 
vortex burst and flow field formation/reattachment 
result in the hysteresis loops seen in the data. These 
lags are seen to significantly delay the unstable 
characteristics both in roll and yaw to higher angles of 
attack than for static conditions. The effect of the 
vertical tail on the dynamic increments in sideslip is 
shown in figure 10. These data show that the vertical 
tail had very little impact on rolling moment 
increments. The unsteady effects did impact the flow 
around the tail giving directional stability 
improvements dynamically compared to the static 
conditions at a = 36°. 

Roll and Yaw Axes Results 

Roll and yaw forced oscillation tests were conducted in 
the NASA Langley 14 -by -22 Foot Subsonic Tunnel 
using the standard low-speed forced oscillation test rig 
previously used in the NASA Langley 30 -by -60 Foot 
Tunnel. Figure 1 1 shows a typical time history of 
rolling moment and roll angle of the F-16XL model 
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oscillation at a = 34° at an amplitude of 30°. Figure 12 
shows rolling moment versus roll angle for 40 
oscillation cycles illustrating the repeatability of the 
dynamic data. Several methods have been proposed and 
used in analysis of forced oscillation data. An 
integrated technique 25 which involves using in-phase 
and out-of-phase components to determine damping 
coefficients is the predominately used technique. 
Another technique, sometimes referred to as the “single 
point method”, utilizes only the points through the 
midpoint of the oscillation cycle, where a and (3 are the 
known nominal values, angular rate is a maximum 
(absolute) value, and angular acceleration is zero. A 
comparison of the two methods was made and is shown 
in figure 13. The single-point values were obtained by 
subtracting the rolling moment coefficient at negative 
rate from the positive rate value, and dividing by the 
sum of the absolute values of the rates. The data were 
obtained at three amplitudes, and the frequency of the 
oscillation was varied to maintain a constant non- 
dimensional roll rate through the midpoint of the 
oscillation cycle of p = 0.04. The results show a large 
effect of amplitude, but both data reduction methods 
resulted in the same values of roll damping. 


the motions, time history effects could be quite 
different. Several other configurations were tested 
similarly to the F-16XL to investigate this possibility. 
Rolling moment data for a modern fighter (fig. 2) with 
low-swept wings, a highly swept LEX, and a long 
forebody are shown in figure 15. This configuration 
shows much less sensitivity to roll oscillation amplitude 
than the F-16XL configuration. Additionally, the trend 
of the data with amplitude is opposite that of the F- 
16XL. Finally, the effects of amplitude persist to the 
extreme a’s instead of decreasing, as was the case for 
the F-16XL. Flow visualization or unsteady pressure 
measurements would be needed to explain the physics 
behind the differences. 

In order to identify key parametric effects, make 
comparisons to current data models, and allow 
comparisons to other techniques, a simple model was 
initially chosen to represent the forced oscillation data. 
The same regression technique discussed previously in 
the pitch results section was used to identify a model of 
the form: 

Cj (a, p, k) = Cj (a,p, k) * pb/2V + C]p(a,p, k) * p. 


Rolling moment data for the F-16XL model at various 
angles of attack and a constant p through the midpoint 
of the oscillation cycle of p = 0.04 is shown in figure 
14. With constant p max values, the points at the 
midpoint of the oscillation cycle have identical state 
values of a, p, and p, and all have zero roll 
acceleration. For angles of attack below 20°, the effect 
of the amplitude of the oscillation is very small, 
evidenced by the coincidence of all of the data at <p = 

0°. The positive values of Cj at cp = 0° correspond to 
negative roll rate, and correspondingly, the negative 
values of Cj at (p = 0° is for positive roll rate. The 

rolling moment increment at (p = 0° between the upper 
and lower segments of each curve is an indication of the 
traditional rate damping. Because tp and P are inversely 
correlated, the slope of the curve is an indication of the 
dihedral effect. 

For angles of attack between 25° and 40°, large 
variation of the data at cp = 0° are seen between the 
oscillations at different amplitudes. This effect is 
analogous to that described previously in pitching 
motions as “time history" effects. Finally, as a is 
increased further, the effect of amplitude is decreased 
and becomes unimportant for this configuration. 

These effects may be very configuration dependent. 
Depending on dominant flow structure topology during 


This model allowed correlation of parameters with k, 
and allowed for the calculation of a parameter to 
compare with conventional rotary balance data. Rotary 
balance data is obtained by rotating the model at 
various angles of attack and sideslip about the wind 
velocity vector. 


The regression model resulted in very good fits of the 
data below a = 35° and above a = 45°. At the 
intermediate angles of attack, the fits were still 
qualitatively good, but higher order terms would have 
to be included in the model to follow the nonlinear 
shape of the data. An example of the regression model 
results is shown in figure 16 for the F-16XL at a = 30° 


and oscillation amplitudes of 5° - 30°. The data show 
that the damping parameter rapidly decreases as k 
increases, while the “static” term (C^) increases with 

increasing k values. Also plotted is the static wind 
tunnel value obtained for reported in a previous 

test 26 . The trend of obtained dynamically appears 

to correlate well with the obtained previously in a 

static test. 


Figure 1 7 shows the regression results for roll 
oscillations at a = 40°. As previously mentioned, this 
angle of attack range produced poorer fits to the wind 
tunnel data, and consequently, more data scatter is 
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observed. Least-squares curve fits to the regression 
model results are shown and indicate the same trends 
seen at a = 30°. Again, the static data previously 
obtained is plotted at k = 0 for reference and the 
dynamically derived values appear to have good 
agreement. Unfortunately, amplitude effects are still 
present, particularly in the mid-angle-of-attack range as 
discussed above. Figure 18 shows the effect of 
amplitude on rolling moment coefficients at a = 40°. 
Amplitude effects are present for yaw oscillation data at 
the higher amplitudes, shown in figure 19. 

Figure 20 shows a summary of the effects of k on 
rolling moment derivatives over a range of angles of 
attack. The data show that at low angles of attack, the 
effects of k are small. At the moderate angles of attack, 
k effects are very large for both Cj and Cjp. As angle 

of attack is further increased, the effect of k decreases 
again. This character of the data indicates that a 
reasonable model may be achievable by including the k 
effects over a fairly narrow range of angle of attack. 

Combining the rolling moments due to roll rate, and 
yaw rate with the formula 

C lQ = C^cos (a) - C[ f sin (a) 

enables a comparison with rotary balance data which is 
obtained by steady rotation (k=0) about the wind 
velocity vector. Since is a steady value, and the 

Cj and Cj^ values are unsteady measurements 

(functions of k), vectorally combining these values, as 
is currently done in some simulation models using the 
Kalviste 4 or other techniques may be misleading if large 
dependencies on k exist for the body axis parameters. 

A comparison of forced oscillation and rotary data is 
shown in figure 2 1 . The rotary data is plotted at k=0, 
and is shown as two boxes on the plot, which covers the 
range slopes of the nonlinear data. 

The intent of inclusion of the Strouhal number is to 
include the effects of frequency content in the data. 

This is expected to account for much of the unsteady 
effects seen in the measured data, and does a reasonable 
job, based on the model fits. Further improvement may 
be realized by analyzing the data retaining 
nonlinearities with roll rate and k. Figure 22 shows an 
example of the rolling moment at constant k values for 
a range of roll rates. These data are obtained by 
oscillating the model at various amplitudes to develop a 
range of rates through the midpoint of the cycles at 
constant k values. The effect of k, and nonlinearities 
associated with rate are seen in this data format. It is 


important to realize that data obtained for this format of 
analysis needs to be at constant k values. The result of 
simply varying oscillation frequency at a constant 
oscillation amplitude results in what could be very 
misleading data, as shown in figure 23, if k effects are 
significant. Using nonlinear forced oscillation data, 
similar to that presented in figure 22, allows for more 
thorough comparisons with rotary data. Although data 
taken with the F-16XL and modern fighter 
configuration presented herein was insufficient for 
these comparisons, data obtained for another 
configuration included sufficient forced oscillation 
results to make comparisons. The results with this 
other configuration showed very similar trends with 
rotary balance data to that seen with the F-16XL 
oscillatory and ramp motion comparison presented 
earlier. At low k values, the trends of the forced 
oscillation derived data agreed well with the rotary 
balance data, and at high k values, large differences 
between the two data sets were apparent. It is 
important to note that, as before, these differences were 
primarily only important in the middle angle of attack 
range around stall. 

Another important effect that has traditionally been 
neglected is the effect of sideslip on body axis damping 
or unsteady terms. Figure 24 shows an example of the 
effect of a non-zero mean angle of sideslip on the roll 
damping parameter (using the previously mentioned 
regression equation) as a function of k at a = 30°. 

These data were obtained by oscillating the model with 
various combinations of sting roll angle offset and pitch 
angles. All data were obtained with oscillation 
amplitude of 30°. The data show that small sideslip 
values have no apparent effect on the coefficient, 
however at larger values of sideslip, large decreases in 
the damping parameter result, eventually even resulting 
in unstable roll rate damping. Additionally, at least 
through the moderate angles of sideslip, large variations 
of the damping parameter with k are seen - similar to 
the zero sideslip data shown previously. 

These results highlight potential sources of errors in 
simulation modeling if the similitude requirements are 
not followed. For example, assume a flight time history 
motion of a periodic rolling maneuver at a = 30° and 
30° maximum roll angle, and k = 0. 1 . Using the 
standard linear superposition model previously 
discussed, and the data in figure 1 6, at the maximum 
roll angle, an error in rolling moment of over 38 % 
(compared to measured wind tunnel data) is produced 
due to not including the effect of k on the “static” term, 
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A limited evaluation of a simple model of roll damping 
- including frequency effects - was conducted as a first 
step in developing more comprehensive and robust 
modeling approaches. The roll damping parameter 
(Cj ) and an increment to the conventional static lateral 


stability (AC [ ^ j was modeled as a function of angle of 

attack and reduced frequency based on the values 
shown in figure 25. The reduced frequency was 
computed by curve fitting a sinusoidal response to the 
time history of roll rate at each time step in the 
simulation. Based on this approach, the frequency was 
computed as; 


w 2 r zMiziil 

~ (P2+Pi)(h ~h) 

and an average value of k for the previous one second 
was used for the table lookup. It was recognized that 
the calculation of frequency from dynamic aircraft 
motion time histories could provide misleading results 
and must be carefully interpreted. However, for the 
purposes of this study, this approach was satisfactory 
for initial evaluation of the aerodynamic effects 
discussed in this paper. Other studies are currently in 
progress that are focused on approaches for modeling 
time-dependent aerodynamic effects and 
implementation in simulation 17,23 . 


Figure 26 shows time history outputs from a real-time 
simulation of a modem fighter airplane using the 
frequency-dependent roll damping model compared to a 
baseline model without the frequency-dependent 
effects. The computed values of k varied from 0 to 0.4 
with a mean value of approximately 0. 1 . The 
simulation results show differences in the motions as a 
result of using the frequency-dependent model 
indicating the importance of the frequency component 
on roll damping, even for relatively benign aircraft 
maneuvers. In addition, this relatively simple modeling 
approach produced smooth and realistic dynamic 
response. Future research is planned to evaluate the 
effects of this modeling approach on flight fidelity at 
high-angles-of-attack. 


Summary and Recommendations 

A consensus appears to be forming that the current 
wind tunnel testing and modeling process is not 
adequate for accurate prediction of flight dynamics of 
modern airplane designs operating at high angles of 
attack and with large angular rates. Many researchers 
are currently developing methods to try and improve on 
the current state-of-the-art. This paper has presented a 
few examples of recent results of large amplitude 
testing of fighter model configurations. The review of 
similitude requirements pointed out the need to 


maintain Strouhal number (reduced frequency) when 
comparing data sets. The effect of reduced frequency 
was configuration dependent, as well as functions of 
other state variables. Direct application of wind tunnel 
oscillatory data will require a method for determining 
an “equivalent” reduced frequency at any given point in 
time to allow the correct data to be looked up. An 
example of a simple representation of reduced 
frequency in simulation time history, and the 
application of a simple rolling moment model using this 
method was shown to be feasible. More complete 
models, including nonlinear effects of rate, and control 
surface deflections, and other state variables may 
require more sophisticated data table look-up 
approaches such as fuzzy logic or neural network 
schemes. 

In order to develop better approaches to modeling the 
unsteady aerodynamic effects, flight validation of the 
models will be essential. This could effectively be 
accomplished using small, unpowered models that 
could be rapidly maneuvered in flight, and with 
sufficient data recording instrumentation to reduce the 
data. These same (or similar) models could then be 
tested using the wind tunnel test techniques available to 
generate the math models and compare directly with 
flight. Once this step is completed, other issues such as 
Mach number, Reynolds number and scale effects 
would need to be addressed. Another important area 
requiring more research, is the methodology of 
combining the various types of data - static, forced 
oscillation, rotary balance, and others which may be 
determined to be required - to approximate the actual, 
non-body axis maneuvers occurring in flight. 
Conventional methods of combining forced oscillation 
and rotary balance data need to be addressed in light of 
the pronounced effects of reduced frequency the forced 
oscillation data may contain. 
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Figure 1 - Photograph of 0.10-scale F-16XL mounted 
for large amplitude pitching tests. 
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Figure 2 - Figure of modern fighter model mounted for 
roll oscillations in the 14x22 Subsonic Tunnel. 
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Figure 3 - Similitude requirements 



Figure 4 - Pitch oscillation results. F-16XL, q max = 0.02, otg = 36°. 



Figure 5 - Effect of pitch rate from ramp motions. F-16XL. 
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qc/2V qc/2V 

Figure 6 - Oscillation amplitude effects. F-16XL, ocq = 21°. 



qc/2V qc/2V 

Figure 7 - Oscillation amplitude effects. F-16XL, ocq = 51° 



Figure 8 - Comparison of ramp and linearized damping parameter, ocq = 36°, F-16XL 
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Figure 1 1 - Time history of rolling moment during roll 
oscillation test. F-16XL, a = 34°. 
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Figure 1 2 - Rolling moment coefficient response with 
roll angle. F-16XL, a = 34°, p max = 0.04. 
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Figure 13 - Comparison of integrated and single point 
methods of computing roll damping parameter. F- 
16XL. 



-32-28-24-20-16-12-8 -4 0 4 8 12 16 20 24 28 32 -32-28-24-20-16-12-8 -4 0 4 8 12 16 20 24 28 32 


<t>, deg <|>, deg 

a =15° a =20° 

a. Low angle of attack 

Figure 14 - Rolling moment response. F-16XL, p max = 0.04. 
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Figure 15 - Rolling moment response. Modem fighter configuration model, p max = 0.04. 
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Figure 1 5 - concluded. 
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Figure 1 6 - Effect of reduced frequency on rolling 
moment derivatives. F-16XL, a = 30°. 



c Ir 


- 0.007 


Figure 1 7 - Effect of reduced frequency on rolling 
moment derivatives. F-16XL, a = 40°. 




Figure 18 - Amplitude effects on rolling moment derivatives during roll oscillations. F-16XL, a = 40°. 



Figure 1 9 - Amplitude effects on rolling moment 
derivative during yaw oscillations. F-16XL, a = 40°. 
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Figure 20 - Effect of angle of attack and reduced frequency on rolling moment derivatives. F-16XL. 



Figure 21 - Comparison of forced oscillation and rotary b. Yaw axis oscillation 

balance results. F-16XL, a = 30. Figure 22 - concluded. 
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a. Roll axis oscillation Figure 23 - Rolling moment results using constant 

Figure 22 - Non-linear rolling moment results. F-16XL amplitude oscillations. F-16XL. 
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Figure 24 - Effect of sideslip on roll damping 
derivative. F-16XL, a = 30°. 




Figure 25 - Roll damping model used for simulation. Modern fighter configuration model. 





Figure 26 - Simulation time-history comparisons. 
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